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Abstract 

^ . We have studied nuclear quadrupole deformation effects in deeply bound 

pionic atoms theoretically. We have evaluated the level shifts and widths 

of the hyperfine components using the first order perturbation theory and 

compared them with the effects of neutron skin. We conclude that the 

. nuclear quadrupole deformation effects for deeply bound Is and 2p states 

are very difficult to observe and that the effects could be observed for 3d 

^ . states. We also conclude that the deformation effects are sensitive to the 

'nT , 

O ■ parameters of the pion-nucleus optical potential. 

O 
l> 

O . PACS: 33.25.Fs, 36.10.Gv 

, 

. Keywords: nuclear quadrupole effects, deeply bound pionic atom 



C ' 1. Introduction 

> : 

^ ■ Deeply bound pionic atoms such as Is and 2p states in heavy nuclei were predicted 

. to be quasi-stable |jl|] and were expected to be observed by proper experimental meth- 

ods 0, Q . There have been a number of efforts to find these states both experimentally 
and theoretically [4-10]. Very recently, deeply bound pionic atoms have been discov- 
ered by Yamazaki et al. [11| using (d,'^He) reactions 12]. In the data we clearly 



identify a peak structure in the energy spectrum of the ^He which is mainly due to the 
formation of deeply bound 2p atomic state with a ^p'^J^ neutron hole state. Formation 
of the deepest Is atomic state is also studied theoretically |T^ and will be tried to find 



experimentally in near future |T3. We can observe the energy levels and widths of 



the atomic states by knowing the reaction Q-value of the peak position and the peak 
width of the ^He spectrum. Hence, we can expect to develop precise spectroscopic 
studies of deeply bound pionic atoms using the reaction with high energy-resolution. 
Since deeply bound pionic atoms exist very close to nuclear surface p|, ||, and since 
they can be populated in all heavy nuclei including [3 unstable ones using the (d,^He) 
reactions in principle 0, we can expect to obtain new information on the nuclear sur- 



face and pion-nucleus interaction from the spectroscopic studies of the deeply bound 
pionic states. 

Nuclear quadrupole deformation effects on pionic processes have been studied 
both in bound states |15|, ^ and scattering states . In the previous works for the 
shallow bound states |]TB|, 0, the quadrupole effects were studied in order to extract 
information on the nuclear deformation using experimental pionic atom spectra. In 
this paper, we study the quadrupole effects for deeply bound pionic states from dif- 
ferent point of view. In order to know the property of the deeply bound pionic atoms 
well and to consider how we can obtain new information from the pionic atoms, we 
need to evaluate the nuclear deformation effects. If the effects are large enough, we 
may be able to use them to observe nuclear deformation precisely, and if the effects are 
small enough, we can just forget them and extract other information like pion-nucleus 
interaction and/or neutron skin. Since the deeply bound states have been observed 
clearly only once, our knowledge is very restricted and we need to study their proper- 
ties in order to develop the spectroscopic study of the deeply bound pionic states. So 
far the deformation effects were evaluated only for the shallow pionic atoms such as 
3d and 4/ in heavy nuclei and were known to cause level splitting and width change 
of the atomic states [1^, |16[. We can expect to have larger deformation effects for 



deeper bound states since pion exists closer to nucleus and, thus, feels nuclear surface 
more sensitively. 

In this paper we evaluate the nuclear quadrupole deformation effects on the deeply 
bound pionic atoms to know how important the effects are in pionic atom spec- 
troscopy. We use the first order perturbation theory to evaluate the level shifts and 
widths for the individual hyperfine components, which is described in section 2. In 
section 3 we show the numerical results. We summarize this paper in section 4. 



2. Quadrupole effects on pionic atoms 



In this section we describe briefly our theoretical formalism which is used to 
evaluate the quadrupole deformation effects on deeply bound pionic atoms. We apply 
the first order perturbation theory |jl5|. 

The nuclear density is defined by the expectation value of the nuclear ground state 
with spin J and its z-projection Jz = J as, 

p(r) =< J J I p(r) I J J >, (1) 

where /3(r) is defined as, 

P(r)=E^(r-r,,), (2) 

1=1 

with nuclear mass number A. By multipole expansion, we can rewrite the density as. 



P(r)=Po(r)+ E \I^^Pdr)Y,\r), (3) 



even k>2 

where for A; > 2 



Here, the k = 2 term corresponds to the nuclear quadrupole density which we will 
consider in this paper. As a practical parametrization, we use the Woods-Saxon form 
for the nuclear density distribution in the intrinsic frame: 

1 + exp(^^) 

with 

R(e)^R(l + (3Y,'(e)), (6) 

where pjv is the normalization constant, /3 a quadrupole deformation parameter. The 
radius parameter is determined as = 1.2 A^^^ fm, and the diffuseness parameter is 
fixed to a = 0.5 fm. For the small deformation, we expand the density in powers of 
/? and take first two terms as. 



p{r)^p-o{r) + J^p-,{r)Y,'{e), (7) 



where 



Pn 

1 + exp(^) 

and 



Poir) - , (8) 



h6^ R exp(^) 
P.{r) = V^/^^^- (i + exp(r^))^ ■ 

The nuclear densities in the laboratory frame p{r) and P2ir), which should be used 
for our calculations for pionic atoms, are related to those in the intrinsic frame as 
follows. 

Pair) = po{r), 



and 



The (3 parameter is related to the nuclear spectroscopic quadrupole moment Q through, 



Q = / pliry'^dr, (11) 

where a superscript 'p' indicates the proton density distribution. We can make similar 
approximation to the p^(r) terms which appear in pion- nucleus optical potential as, 



p'(r) 



pl{r) + 2J—p,{T)p,{T)Y^{e). (12) 



We can calculate pion wave functions by solving the Klein-Gordon equation with 
the Coulomb potential of a finite nuclear size, and with the pion-nucleus optical 
potential. We include only spherical nuclear densities for the Klein-Gordon equation 
and evaluate the quadrupole deformation effect by perturbation. We have applied 



the standard Ericson-Ericson potential with a parameter set determined by Seki 



et al. [0 and Konijn et al. [^. All the details of the calculation are described in 
Ref. i. 

Using the (3 parameter, the potentials can be written as 

Vopt — ^opt "T H <'opt ) 

and 

Vc = V^'^ + PV^'\ (13) 

The second terms are due to quadrupole deformation and treated as perturbation. 
Then, we can express the complex energy shift in the following way: 

AE = p[<vSI> + <^^vi'^>], (14) 

p 

where E is the complex eigenvalue of the Klein-Gordon equation and p the reduced 
mass of pion. The terms coming from Coulomb and optical potential contribute to 
the complex energy shift differently as shown in Eq. (14) since these two potential 
appear in the Klein-Gordon equation in a different way as in Eq. (17). 

Because the Klein-Gordon equation includes the complex optical potential which 
makes the Hamiltonian non-Hermite and makes the eigen energies complex, we should 
be careful to the normalization condition to get proper orthonormal wave functions. 
The Klein-Gordon equation can be rewritten as the following two-component equation 



21 



((T3 + tcr^)^ + a^p + (as + 2^2)^?* + V^c = E^-^i, (15) 



where (Tj's are the Pauli spin matrices and 



This equation is equavalent to the Klein-Gordon equation, 

[-V^ + /i^ + = (E, - V^y<P,. (17) 

/^From Eq. (15), we get orthonormal condition for \1/ as; 

^fas^jdr = 5ij, (18) 



where the \l/ means to take complex conjugate only for the angular part |21|]. This 
provides proper orthonormal condition for the wave functions as follows; 



0. 



(0) 

dr = Sij (19) 



where the means to take complex conjugate only for the angular part of the wave 
function. We have used the wave function, which satisfies Eqs.(17) and (19), to 
calculate all expectation values in this paper. 

Then, we evaluate the deformation effects on pionic atoms using the first order 
perturbation theory. To calculate the complex energy shift in Eq. (14), we need 
following expectation values for the deformation effect due to strong interaction: 

^ 7 i-oo 

< p(r) >U=2= -^C{JIF)-— / Rr,,{r)p,{r)RUrydr, (20) 
87r 21 + 6 Jo 

+ [1(1 + I) -3] Rni{r)p2{r)Rni{r)dr}, (21) 
Jo 

5 ^, / 



< p\v) >U=2= -^^'(J/F)^ RUr){2po{r)p,{r))R^,{rydr, (22) 



and 



<Vp^(r)V + ^C{JlF)-J—{ r 

Sn 2/ + 3 JO dr dr 

+ [/(/ + l)-3]/ Rniir){2po{r)p2ir))Rni{r)dr}. (23) 



For the deformation effects coming from Coulomb interaction, we need the matrix 
element: 



k=2 



roc POO /r-^ 

r-r'l — 2 2/ + 3/0 '"''^'"/o ^"d^' ^Rni{r)Rni{r)p'2{r'). 

(24) 

In above expressions, Rni{r) is the radial wave function of the bound pion and we 
have used the coefficient C{JIF) defined as, 

^ 3X(X-1)-4J(J+1W + 1) 

^ ^ 2J(2J- l)/(2/- 1) ' ^ ' 

with 



X = J(J+l) + /(/ + l)-F(F + l), (26) 

where / is the angular momentum of pionic states, and F = I ® J the total angular 
momentum of the system. Obviously, the quadrupole deformation has no effect on 
pionic s states within the first order perturbation. 



3. Numerical results 



In this section, we evaluate the nuclear quadrupole effects on the deeply bound 
pionic atoms. We show the results obtained with the potential parameter set [a-1] 
in Table 1 of ref. In the numerical calculations, we have used the computer 



program MATOM |]T9[ to obtain eigen energies and wave functions of pionic 2p and 
3(i states. As a typical example of the deformed nucleus, we consider ypLuio4 whose 



ground state is J'" = ^ and quadrupole moment Q = +3.49 (b). We fix the density 
parameters appeared in the previous section as i? = 6.71 fm and a = 0.5 fm which 
give jS = 0.29. By solving the Klein-Gordon equation with the spherical density, 
we obtain the binding energies and widths of the pionic 2p and 3d states which are 
B.E.(2p)=4084, r{2p) = 258 , B.E.(3rf)=2133 and r(3d) = 29 in unit of keV. 

First, we evaluate the quadrupole deformation effects on pionic 2p and 3d states. 
In Fig. 1, we show the level shifts and widths of the pionic 2p state (a) and 3d state (b) 
as a function of the (3 parameter. We can see in the Fig. 1 (a) that the 2p level splits 
into three levels with spacing of order lOkeV, which is one order of magnitude smaller 
than the width of the 2p state due to the nuclear absorption of bound pion. Thus 
the quadrupole hyperfine splitting will not change the level structure significantly 
for this state. The 3d state, as shown in Fig. 2 (b), splits into five levels whose 
spacing is the same order of magnitude as the width of the 3d state. In this case, 
the level splitting could be seen by proper experiments. For example, we could see 
the levels as small separated peaks on the top of the large peak consist of all 3d state 
contributions in double differential cross section of the (d,^He) reactions, if we can 



achieve high energy-resolution hke 10 keV. Even with worse energy resolution, we 
may be able to see the effects by observing the significantly larger width compared 
to those calculated with spherical nuclear density. 

Next, we evaluate the quadrupole effect for the nucleus with neutron skin. This 
evaluation is much interesting since deeply bound pionic atoms are known to be 
sensitive to the neutron skin ||2^ and since the unstable nuclei are expected to have 
the neutron skin and the deformation simultaneously. To simulate the existence of 
the neutron skin, we introduce a parameter A as. 



A = -R„ — Rp, (27) 

and vary the A by changing the radius of neutron distribution Rn. We show the 
calculated energies of deeply bound states as a function of A in Fig. 2. As shown in 
the figure, the deepest Is state is most sensitive to the thickness of the neutron skin. 
We also show the calculated energy levels of 2p and 3d atomic state as a function of 
deformation parameter /3, in Fig. 3 for A = 0.5 fm and in Fig. 4 for A = 1.0 fm. 
These figures show similar behavior as in Fig. 1 and indicate that the deformation 
effects on pionic states are qualitatively same for nuclei with thick neutron skin. 

It is also interesting to compare the magnitude of the deformation effects with 
neutron skin effects. For this purpose, we consider two cases A(fm)] = [0.5, 0.0] 
and [0.0, 1.0], and compare the calculated energies in each atomic state. The results 
are foUowings: 

The deepest Is state is very sensitive to the neutron skin. The energy shift and 
width change due to neutron skin are \AE\ ~ 300 keV and |Ar| ~ 900 keV, while 
the energy of the Is state is not affected by the deformation within the first order 
perturbation theory as mentioned at the end of section 2. 

The sensitivities of the 2p state to neutron skin are \AE\ ~ 200 keV, |Ar| ~ 10 
keV. Those to deformation are \AE\ = 10 ~ 40 keV, |Ar| = 10 ~ 70 keV. The 
deformation effects depend on the total spin F of the system defined in the last 
section. 

The sensitivities of the 3d state to neutron skin are \AE\ ~ 20 keV, |Ar| ~ 10 
keV. Those to deformation are |AE| = ~ 20 keV, |Ar| < 10 keV. 

Before closing this section, we describe numerical results obtained with another 
set of potential parameter determined by Konijn et al. We used the parameter set 
'present fit, ^ = V shown in Table 1 of ref. |j20|. We have calculated the binding 
energies and widths of 2p and 3d state with the spherical nuclear density and obtained 
the similar results as those with Seki's potential |]19|. At A = [fm], we obtained 
B.E.(2p) = 4107, r(2p) = 257 and B.E.(3d) = 2140, T{3d) = 31 in unit of keV. 

Then we evaluated the quadrupole deformation effects with Konijn's potential and 
compared them with those by Seki's [a-1] potential. For 3d state, Konijn's potential 
gives 40-60% larger quadrupole shifts with the same sign. As for the width, we 
obtained much larger deformation effects and its signs are opposite to those with 



Seki's potential at A = and 0.5 fm, while both potentials give the quadrupole 
widths with the same sign for A = 1 fm. For 2p state, Konijn's potential gives 
opposite sign for both the real and the imaginary parts of the quadrupole energy 
shifts. Concerning the absolute values, we obtain larger real energy shifts and smaller 
imaginary energy shifts with Konijn's potential for A = and 0.5 fm, while for A = 1 
fm Seki's potential yields larger real energy shifts and smaller imaginary energy shifts. 

We have found that the deformation effects are sensitive to the potential parame- 
ters because of the following reasons: As was pointed out by Seki and Masutani jTotl , 
the potential parameters are correlated with each other. The isoscalar s-wave term is 
approximately written as (bo + Bope)p where pe is the effective nuclear density. Then 
the potential parameters with different bo and Bo but with similar (bo + BoPe) value 
yield almost the same level shifts and widths for the pionic atom for spherical nuclei. 
This is the same for p-wave terms. In fact, the real part of the potential parameters 
by Seki [a-1] and Konijn = 1] are quite differnt but the numerical values (bo + Bope) 
are almost the same. For the deformation effects, however, the different combination 
of the potential parameter appear as < (6o + 2-Bop)p2 >~ (&o + 2-BoPe) < P2 > and 
these are considerably different for the potentials by Seki [a-1] and Konijn. Thus, 
the potentials which yield the similar results for the level shifts and widths of the 
spherical nuclei possibly give different quadrupole contribution. Futhermore, there 
are large cancellation between the deformation effects coming from the strong and 
Coulomb interaction, so that the moderate changes of the strong quadrupole effects 
gives the large modification to the total deformation effects. Also, because of the 
complex pion wave function, the Coulomb term has contributions to real and imag- 
inary energy shifts. These are the reasons why the potentials by Seki and Konijn 
predicts fairly different quadrupole effects for deeply bound pionic atoms. These situ- 
ation does not happen for the quadrupole effects in shallow pionic states, because the 
Coulomb contribution dominates and there are no large cancellation between strong 
and Coulomb contributions. 

4. Summary 

We have evaluated the nuclear quadrupole deformation effects on deeply bound 
pionic atoms. The nuclear deformation induces the hyperfine splitting and affects 
the widths of the pionic states, which are expected to be larger for deeper states 
since pion exists closer to the nucleus. We are very much interested in knowing the 
properties of the deeply bound pionic atoms for deformed nuclei and to find out the 
best atomic state - nucleus combinations for obtaining new information clearly. 

We found that the level shifts and width changes due to the quadrupole defor- 
mation are too small to observe for deeply bound Is and 2p states. For 2p state the 
natural width of the state due to pion absorption is one order of magnitude larger 
than the level shift due to nuclear deformation. There are no quadrupole effects on 



pionic Is states in the first order perturbation. Thus, it is very unhkely to observe the 
quadrupole effect in pionic Is and 2p states. For atomic 3d state, we have a chance to 
see the level splitting due to the deformation since the natural width of the 3d state 
is only around 30 keV, which is comparable to the level shift due to the quadrupole 
deformation. Thus, they may be observed with proper experimental technique like 
(d,^He) reaction with high energy-resolution. 

We have also studied the quadrupole effects for nuclei with neutron skin. This is 
interesting since we expect to have skin and deformation simultaneously in unstable 
nuclei, where we can produce deeply bound pionic states in principle using the (d,^He) 
reaction in inverse kinematics p. We change the radius of the neutron distribution 
artificially to simulate the neutron skin. It is found that the existence of the neutron 
skin does not change the deformation effect significantly. We compare the magnitudes 
of the both effects and summarize the results in section 3. 

We have evaluated the deformation effects using two different optical potential 
parameters and found that the nuclear quadrupole deformation effects are very sensi- 
tive to the potential parameters in deeply bound pionic states. Thus, we can expect 
to determine the optical potential parameters more precisely if we can observe the 
nuclear deformation effects on deeply bound pionic atoms systematically. 

/^From the results in this paper, we can conclude that the nuclear quadrupole de- 
formation effects are difficult to observe experimentally for deeply bound pionic states 
like Is and 2p orbits. In order to determine the deformation from the pionic atoms, 
which could be interesting for unstable nuclei, we need to determine neutron skin 
thickness (monopole nuclear density) firstly by observing the Is state and then de- 
termine the deformation from the observation of 3d state with high energy-resolution 
like 10 keV. 
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Figure Caption 

Fig. 1 Binding energies of (a) pionic 2p state and (b) 2>d state as a function of the 
deformation parameter /3 defined in the text. Total spin F of each state is indicated 
in the figures. The width of each state is shown by the error bars. 

Fig. 2 Binding energies of pionic Is, 2p and 3(i states as a function of neutron skin 
thickness A defined in the text. The width of each state is shown as hatched areas. 
The hatched area of 3d state is obtained by multiplying a factor 10 to the calculated 
width. 

Fig. 3 Same as in Fig. 1 except for neutron thickness parameter A = 0.5 fm. 
Fig. 4 Same as in Fig. 1 except for neutron thickness parameter A = 1.0 fm. 
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